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Abstract 

In an interesting recent paper on the growth of inhomogeneity 
through the effect of gravity [1], Bertschinger and Hamilton derive 
equations for the electric and magnetic parts of the Weyl tensor for 
cold dust for both General Relativity and Newtonian theory. Their 
conclusion is that both in General Relativity and in Newtonian the- 
ory, in general the magnetic part of the Weyl tensor does not vanish, 
implying that the Lagrangian evolution of the fluid is not local. We 
show here that the 'Newtonian' theory discussed by them is in fact 
not Newtonian theory per se, but rather a plausible relativistic gen- 
eralisation of Newtonian theory. Newtonian cosmology itself is highly 
non-local irrespective of the behaviour of the magnetic part of the 
Weyl tensor; in this respect the Bertschinger-Hamilton generalisation 
is a better theory. 



1 The context 



The question at issue [1] is the relation between the fully covariant analysis 
of fluid flows in general relativity [2-5] as compared with that of Newtonian 
gravitational theory, and their implications for the study of inhomogeneities 
in cosmology ([6-9] and references therein). The particular point that has 
generated interest has been the recent discovery of a family of General Rela- 
tivity exact solutions - those with vanishing magnetic part of the Weyl tensor 
- where neighbouring flow lines evolve independently of each other, in the 
sense that this evolution depends only on the local fluid variables [10-14]. 
Actually it has been known for a long time that this is true for spherically 
symmetric solutions, but the solutions now under consideration include the 
relativistic version of the Zeldovich theory of gravitational collapse [15], how- 
ever predicting filamentary gravitational collapse rather than sheetlike grav- 
itational collapse as in the Zeldovich analysis. It has been suggested that 
this difference arises because Newtonian theory is essentially non-local, and 
that the relativistic theory should have a non-zero magnetic part of the Weyl 
tensor to give corresponding non-local physics. The issue then arises, does 
the magnetic part of the Weyl tensor necessarily vanish in the Newtonian 
limit (as claimed in [4]), and are the Lagrangian evolution equations in that 
limit local in general? According to [1] , the answers to both questions are no. 

However this discussion has not referred back to the analyses of New- 
tonian cosmology by Heckmann and Schucking [16-18], where it is shown 
firstly that no Newtonian cosmology - based on a potential and Poisson 
equation - is possible without some extension of Newtonian theory as nor- 
mally understood; and that the obvious extension of Newtonian theory to 
the cosmological situation is essentially non-local - so that local physics can- 
not be decoupled from instantaneous boundary conditions at infinity. This 
is because Newtonian theory is a singular limit of general relativity theory, 
precluding the possibility of gravitational waves proceeding at a finite speed. 
It does not in an obvious way involve an analogue of the magnetic part of 
the Weyl tensor [4] . 

By contrast, the Bertschinger-Hamilton theory [1] leads to influences 
propagating at a finite speed, and so their effect do not arrive at any event 
instantaneously from infinity. In this sense their theory - which allows grav- 
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itational waves - is a better approximation to the correct classical theory 
of gravity (namely, general relativity) than is Newtonian theory proper (or 
more strictly, the Heckmann-Schucking version of Newtonian theory [16-18], 
with boundary conditions allowing spatially homogeneous cosmologies). In 
what follows, to save repetition of these names we shall refer to the different 
theories as BH (Bertschinger-Hamilton [1]), HS (Heckmann-Schucking [16- 
18]), NG (standard Newtonian gravitational theory expressed in a potential 
formalism), and GR (general relativity theory). 

2 General Relativity Theory 

In GR applied to cosmology in a covariant manner [2-5], one decomposes ge- 
ometrical and physical quantities relative to a preferred 4-velocity vector u"" 
- the fundamental velocity that underlies any coherent cosmological model 
[2-5]. This defines the fluid expansion 9, shear Uab, vorticity LOab-, and ac- 
celeration tta = iia- Their time derivatives (denoted by a dot, e.g. ©) are 
determined by the matter density /i and pressure p (here we restrict our con- 
sideration to perfect fluids) together with the electric part Eab of the Weyl 
tensor, whose derivative in turn is determined by the magnetic part Hab of 
the Weyl tensor, where these quantities are deflned respectively from the 
Weyl tensor Cabcd by 

Eac = CabcduV, Hac = CabefU^V^^ chu'' (l) 

(each being a trace- free symmetric tensor that is orthogonal to u"). The 

quantities Eab and Hab obey equations very similar in structure to Maxwell's 
equations written relative to a general family of observers [3-5] . A particular 
consequence is that Eab and Hab each obey a wave equation, with somewhat 
complicated source terms (sec [3] for the case of an almost-Robertson- Walker 
space time) . 

In the case of interest, we specialize to 'dust', that is the pressure p 
vanishes and consequently the fundamental flow lines are geodesic (oa = 
iia — 0). The resulting equations are given below. 



3 



3 Newtonian Theory and Cosmology 



Attempts to use Newton's force law in the case of an infinitely spread out 
medium, such as in a spatially homogeneous Newtonian cosmology, are plagued 
by infinities (because the matter source of the gravitational force is un- 
bounded), so in setting up the equations of Newtonian cosmology, an ap- 
proach based on a potential $ is preferable. 

The potential must satisfy the Poisson equation: 

= -AixGp (2) 

where p is the matter density. This has to be supplemented by boundary 
conditions at infinity in order to get a unique solution; the usual ones are 

lim $ = . (3) 

r— >oo 

We take equns. (2,3) as defining standard Newtonian gravity. Now the prob- 
lem is that these cannot give us a sensible cosmological model, for (3) is 
incompatible with a homogeneous distribution of matter [16]; in order to 
deal with an infinite matter distribution we have to drop (3) and replace it 
by some other boundary condition which allows the potential to diverge at 
infinity, as for example in the Newtonian version of the Robertson- Walker 
universe models. 

What are suitable boundary conditions? Following HS, we define the 
Newtonian analogue of Eah by the equation 

E^p = ^^^p-hi^pVH E^p = Ep^, E"„ = 0. (4) 

Then the Newtonian version of Robertson- Walker models are compatible 
with the condition 

lim = , (5) 

r— »oo ^ 

which is thus a weakening of condition (3) that allows one to handle spatially 
homogeneous cosmological models, which necessarily have uniform density^. 

^One also has to generalise the idea of inertial motion to free-fall motion, in analogy 
with general relativity, in order to define spatially homogeneous Newtonian models such 
as the Newtonian version of the FRW models; that issue is not of concern here. 
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However this conditions excludes the Newtonian analogs of the Kantowski- 
Sachs and Bianchi spatially homogeneous but anisotropic universe models - 
which are amongst the simplest generalisations of the FRW universes. Thus 
HS proposed instead the condition 

lira = Eai3{t)\oo (6) 

where Eaf^{t)\oo are arbitrary functions of time 0. Thus the limit of the 
components Eajj at spatial infinity may be arbitrarily prescribed as functions 
of time; they then propagate in to any point with infinite speed, according 
to the equation 

E°'i3,a = -^Gp^a (7) 

which follows from the Poisson equation (2) and definition (4). 

There has been some debate over these boundary conditions, Narlikar 
[19] suggesting that (5) are in fact better conditions to use for Newtonian 
cosmology than (6) - but thereby excluding from his consideration that fam- 
ily of NT models that correspond to the Bianchi GR solutions allowed by 
(6) ([17,20]). Both are of course generalisations of the 'true' Newtonian con- 
ditions (3). Incidentally, one can ask here why one does not specify either 
lim^^oo ^ = ^{t)\oo or lim^^oo V^$; the reason is that the first diverges, while 
the second is determined by the limiting mass density at infinity (through 
the Poisson equation (2)). 

The point that is fundamental to us here is that NT itself gives no equa- 
tion for and consequently can give no equation for E'^/s (with Ea/3 defined 
by (4)). This is true whether one writes NT in Lagrangian or Eulerian co- 
ordinates. Any specification one may obtain for these time derivatives thus 
either results from positing a theory that is not NT itself, but rather some 
generalisation of NT; or from positing some set of boundary conditions, such 
as (5) or (6). The latter option does not directly give any local equation for 
E'ap, which is then determined non-locally by the chosen condition at infinity 
plus the local matter conservation equations and the divergence equation (7) 

^in general they should depend on polar angles 9 and </> also; however in the context of 
spatially homogeneous analogues of the Bianchi models, they can depend only on t. 
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for Ead- NT itself (which obeys (3)) gives no clue as to which of these gen- 
eralised boundary conditions ((6) or its specialisation (5)) we should choose. 

We support the HS view that (6) rather than (5) is the better option. In 
both cases, information immediately propagates in from infinity to determine 
the local physical response; the fact that the information imparted in (5) is 
the 'null' case - there is no change in E^js at infinity - does not change the 
fact we are determining what happens locally by a choice of conditions at 
infinity; and according to (6), these are arbitrarily specifiable. 

One might ask what are the suitable conditions to use at infinity for a 
Newtonian version of a perturbed FRW model, that has statistically spatially 
homogeneous inhomogeneities superimposed on a FRW background (which 
is necessarily spatially homogeneous and therefore stretches to infinity). It 
would seem we are in trouble here, for none of the conditions above seem 
adequate to this case - because then there will be no regular limit for EaB as 
we go to spatial infinity (because of the statistical fluctuations in the matter, 
this quantity too may be expected to fluctuate statistically). We can get 
a good description either by imposing periodic boundary conditions, as is 
done in most numerical simulations - thus avoiding the problem in the way 
first suggested by Einstein through his proposal of spatially compact universe 
models (and corresponding to the possibility of finding Newtonian versions of 
the 'small universe' idea [21]); or by insisting that the perturbations die away 
outside some bounded domain, thus allowing (5) for example as a boundary 
condition, at the expense of denying the assumption of spatial homogeneity 
of conditions in the universe - which is one of the central tenets of current 
cosmological dogma. 

Any use of NT to discuss the growth of perturbations in a cosmological 
context should make it quite clear which of these options is adopted - and why. 
It may perhaps be claimed that physically these both give reasonable results 
- but that has to be shown. In our view the best alternative would be a third 
option that has not so far been systematically developed: namely, to adopt 
a Newtonian version of the 'Finite Infinity' proposal for the GR case [22], 
that is, to isolate the considered local system by a sphere that is far enough 
away to be regarded as infinity for all practical purposes but, because it is 
at a finite distance, can be investigated easily and used as a surface where 
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boundary conditions can be imposed (and the residual influence of the outer 
regions on the effectively 'isolated' interior can thus be determined). 

4 The Bertschinger-Hamilton Theory 

The BH theory is developed as follows [1] : the continuity and Poisson equa- 
tions are written in 'comoving' coordinates, with a background FRW expan- 
sion factored out; thus these equations are written for the fractional density 
perturbation 6 and an associated part (p of the full potential $. This hides 
a tricky gauge problem, because 'comoving' here actually means relative to 
a fictitious background FRW space-time (the matter moves relative to the 
chosen frame) ; one might suggest it would in this context be preferable to use 
the Newtonian version [23] of the GR gauge-invariant and covariant theory 
[6-9]. In the BH approach, the Newtonian gravity vector g is defined from 
the perturbative potential (phj g = — V0, and then in turn defines the tensor 
Eij by a perturbative version of (4), written in comoving coordinates. 

BH then develop a series of equations that are very similar to those ob- 
tained from the linearised GR equations. Substituting g into the continuity 
equation gives 

V. (^^) = 4.Ga3V./, f^pv^f,+h (8) 

where the mass current has been decomposed in the comoving frame into 
longitudinal and transverse parts obeying 

Vx/j| = 0, V./1 = (9) 

— * 

Then the transverse mass current is replaced by a transverse vector field H 
obeying 

V X H = -levrGaVl, V.# = (10) 
From this in turn they define a traceless tensor 

Hij = —VjHi + 2vke''\V^gi = H^, + e^j^A" (11) 
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where Hij is the symmetric part of Tiij and the dual of the anti-symmetric 
part. Now BH perform a series of manipulations, the key ones of which are 
integrating (8) and substituting to get 

^+^^=477^0^+^ (12) 

dr a 

and then taking the spatial gradient of this equation and substituting to get 

—^+'-Eij—Wk^^\iHj)i+QEij—hija^^Eki—'ia^{iEj)k—^^ = —AnGa^paij 

(13) 

At this point there is something of a non-sequitur where a GR equation for Eij 
is introduced; taking the curl of this equation and making some substitutions 
BH obtain the time-derivative equation for Hif 

+ ~^i3 + '^k€^\iEj)i + ^Hij — hija^''Hki — 3(7*^(jijj)fc —uj^{iHj)k = 0- (14) 

Together with the divergence equations for E^j and Hij, these give the set of 
Maxwell-like equations for E^j and Hij that occur in the linearised version of 
GR. 



At this point it is quite clear that we have a theory that is something 
other than NG, because (1) NG cannot determine the time evolution of Eab 
by an equation like (13), as discussed in the previous section, and (2) if we 
substitute (14) into the time derivative of (13) we get a wave equation for 
Eab, showing the possibility of gravitational waves (with speed unity in the 
chosen units) in the BH theory (as in the linearised GR theory, see [3]). How- 
ever it is also clear on comparison with the GR equations (see [1]) that the 
BH theory is a good generahsation of Newtonian theory, in that it gives a 
good approximation to the GR equations. 

How has it come about that the BH theory has produced a time evolu- 
tion equation for Eij, and hence (with the time evolution equation for Hij) 
a wave equation for Eijl The first point is that, in going from (8) to (12), 
the solutions to the equations are very under-determined, for example, there 
is considerable arbitrariness in g. Again in the series of manipulations that 
leads from f to H to Hij and A^, there is considerable arbitrariness in these 
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quantities if one allows for the general possible solutions. In fact if one traces 
what happens in these equations, actually (from the Newtonian viewpoint), 

— * 

the tensor Hij is arbitrarily specifiable. This is because / determines only the 

— * — * 

curl of H (by (10)); on integrating to determine H, the trace-free symmetric 

— * 

part of H is (at least locally) arbitrarily specifiable; but this is just Hij. More 
precisely, the integrability equations for Hi to exist, given Hij, arc just the 
'divergence' equations for Hij ((43) in BH); so Hij can be chosen as an arbi- 
trary solution of these spatial equations, allowing arbitrary time evolution. 
Furthermore because of the arbitrary functions that occur in integrating (8), 
there is also considerable arbitrariness in A\ 

How then do we get an equation ((13) here, (49) in BH) that gives the 
time derivative of Eij in terms of Hijl Our claim is that in fact in this 
equation, both dEij/dt and Hij are arbitrarily specifiable in NT. Prom this 
viewpoint, Eij is arbitrarily specifiable as a function of time, and (13) then 
tells us what Hij has to be, given the definition (11) (or one can run this 
the other way: choose Hij as you like, and (13) then determines dEij/dr). 
Hence this equation does not actually determine the time evolution of Eij - 
for there is no Newtonian equation for the time evolution of Hij (in a truly 
NT approach, we might define H^j by (13), and accept whatever identities 
follow from this definition; none then determine the time evolution of the 
system [17]). 

What then of the BH equations that determine the time derivative of 
Hij (and hence lead to the wave equation already commented on)? Well, in 
order to derive them BH introduce their equation (52), which follows from 
GR rather than NT. Thus it is here that they close the equations in causal 
terms - by importing relations that do not in fact follow from NT. The time 
derivative equation (14) for Hij (their equation (55)) then follows, and hence 
the existence of gravitational waves (and the implied speed of travel of those 
waves) . 

In summary, the Eij evolution equations arise from use of definitions that 
bring the NT theory into the form of the GR equations, but do not in fact 
determine the time evolution uniquely, as long as we remain within the con- 
fines of NG (rather these equations define the variable Hij which can absorb 
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any chosen time evolution). This does not prove wrong the set of equations 
given, which in fact allow rather arbitrary solutions; however the appearance 
of uniqueness is misleading, resulting from non-uniqueness of solutions to the 
equations presented. One can obtain uniqueness by use of special rather than 
general solutions to the equations (just as (5) is a special case of (6)), but 
this is an arbitrary restriction on the allowed solutions of the theory. The 
Hij evolution equation does not result directly from NT at all, but rather is 
imported from GR. 

Thus the BH theory becomes determinate locally (and in good accord 
with linearised GR) only by introducing extra equations which do not in fact 
follow from NG. The result is a theory that is a better approximation to GR 
than NT is, in particular because it allows gravitational waves. This theory 
is more local than Newtonian theory proper, as it has Cauchy development 
properties like GR. It also has the advantage of having a quantity Hij that 
corresponds in the appropriate way (in terms of its role in the equations) to 
the magnetic part of the Weyl tensor in GR (NG theory proper has no such 
tensor [4]). 

5 The better theory? 

Which then of these theories is the better theory? Our viewpoint is that the 
best is GR because it is the most fundamental of all these theories - indeed 
it is the fundamental classical gravitational theory. The true theory is GR. 
The issue is what are usable approximations in Newtonian-like conditions 
(for example, in studying local astrophysics). 

BH is a good candidate, but is certainly a generalisation of NT rather 
than being classical NT. It reflects the causal structure of GR, but is closer 
to hnearised GR rather than full GR (this is no accident; it was constructed 
that way). Prom this viewpoint, NG in turn is an acceptable approximation 
to BH theory in some circumstances; but the range of conditions adequately 
covered by BH theory is wider than that of NG but less than that of GR. 

Which is more useful in astrophysical contexts will depend on those con- 
texts. Thus we have not considered here the issue of which theories give 
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adequate description of anisotropic collapse situations leading to the forma- 
tion of structure in the expanding universe. However the line of argument 
above suggests that where NT and BH theory disagree, we should believe 
the latter (but where BH and GR disagree, we should beheve GR). 

Finally does Newtonian theory itself demand a non-zero magnetic Weyl 
tensor analogue? No - it is probably not even well-defined. But that fact does 
not undermine the claims of BH theory to give a better description of what 
will happen in some collapse scenarios. But that does not mean we should 
necessarily expect a non-zero Magnetic Weyl tensor (or its quasi-Newtonian 
analogue) in realistic collapse situations. The point is that one should treat 
carefully the claim that neighbouring world-lines do not influence each other 
in 'silent universe' [11], for they do indeed feel each other's gravitational in- 
fluence [24] through the set of constraint equations - including the relativistic 
version of the Poisson equation - that are necessarily satisfied in a full solu- 
tion of the field equations (and arc consistent with the evolution equations 
in such 'silent' universes [25]). It is gravitational induction and gravitational 
waves that are precluded in these solutions - but these effects are not likely 
to be important during the Newtonian-like phase of gravitational collapse. 

We thank the FRD (South Africa) for financial support. 
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